Branched Polymers 

Richard Kenyon*and Peter Winkler ' 



Abstract 

Building on and from the work of Brydges and Imbrie, we give an elementary 
calculation of the volume of the space of branched polymers of order n in the plane 
and in 3-space. Our development reveals some more general identities, and allows 
exact random sampling. In particular we show that a random 3-dimensional branched 
polymer of order n has diameter of order ^Jn. 

1 Introduction 

A branched polymer of order n in M. d — or just "polymer" for short — is a connected set of 
n labeled unit spheres with nonoverlapping interiors. We will assume that the sphere labeled 
1 is centered at the origin. See Figured] for an example in the plane. 

Intended as a model in chemistry or biology, branched polymers are often modeled, in 
turn, by lattice animals (trees on a grid); see, e.g., 0, El [7J QUI [TBI EH] . However, we will see 
that continuum polymers turn out to be in some respects more tractable. 

The set of polymers can be parametrized locally by the spherical angles of the vectors 
connecting adjacent sphere centers. In these coordinates Brydges and Imbrie [2] showed 
that the space B d (n) of polymers of order n has total volume (n— l)!(27r)™~ 1 for d = 2 and 
ri n_1 (27r) n ~ 1 for d = 3. Their proof uses nonconstructive techniques such as equivariant 
cohomology and localization. 

We give here an elementary proof, together with some generalizations and an algorithm 
for exact random sampling of polymers. In the planar case our algorithm has the added 
feature of being inductive, in the sense that a uniformly random polymer of order n is 
constructed from one of order n—1. 

Although not explicit in their paper, the proof in [2] in fact shows that in the planar case 
the volume of the configuration space is unchanged when the radii of the individual disks 
are different. We use this fact in an essential way in our constructions. 
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Figure 1: A branched polymer in the plane 

2 The Planar Case 

Let us observe first that (n— l)!(27r) n ~ 1 is also the volume of the space of "crossing worms"— 
that is, strings of labeled touching disks, beginning with disk 1 centered at the origin, but 
now with no constraint that disks may not overlap. See Figure |2] below for an example. 
Fixing the order of disks 2 through n in the crossing worm yields an ordinary unit-step walk 
in the plane of n—1 steps. 

Yet another space of volume (n— l)!(27r) n ~ 1 is the space of "crossing inductive trees", one 
of which is illustrated in Figure [31 A crossing inductive tree is a tree of n touching labeled 
disks with overlapping permitted, but required to satisfy the condition that for each k < n, 
disks 1, . . . , k must also form a tree. In other words, the vertex labels increase from the root 
1. We will see that this space is in fact a certain limiting case of the space of polymers. 

2.1 Coordinates 

In the volume calculation we will need to consider polymers made of disks of arbitrary radius. 
Let Ti € (0, oo) be the radius of the ith disk and R = {r 1; . . . ,r n } be the vector of radii. 
Given a polymer X = X(R), define a graph G(X) with a vertex for each disk of X and an 
edge between vertices whenever the corresponding disks are adjacent. Almost surely G(X) 
is a tree, that is, has no cycles. When G(X) is a tree, we root G(X) at the origin, and 
direct each edge away from the origin. This allows us to assign an "absolute" angle (taken 
counterclockwise relative to the X-axis) to each edge. Let ex, ... , e n -\ be the edges (chosen 
in some order) and 9\, . . . , Q n -\ the corresponding angles. 

For a given combinatorial tree T, the set of polymers X = X(R, T) with graph G(X) = T 
can thus be identified with a subset of [0,27r) n_1 . Call this set BPr(T). The boundary of 
BP^(T) corresponds to polymers having at least one cycle; the corresponding plane graphs 
G(X) are obtained by adding one or more edges to T. Indeed, the boundary of BP R (T) is 
piecewise analytic and the pieces of codimension k correspond to polymers with k (facial) 
cycles. 
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Figure 2: A crossing worm 




A polymer X with cycles lies in the boundary of each BPr(T) for which T is a spanning 
tree of the graph G(X). Each such BP^(T) will contribute its own paramaterization to X. 
Note, however, that some of the spanning trees may be unrealizable by unit disks (e.g. the 
star inside a 6- wheel); we just regard them as BP^(T) of zero volume. 

We can construct a model for the parameter space of all polymers of size n and disk radii 
R by taking a copy of BP/j(T) for each possible combinatorial type of tree, and identifying 
boundaries as above. Note that the identification maps are in general analytic maps on 
the angles: in a polygon with k vertices whose edges have fixed lengths r 1; . . . , r^, any two 
consecutive angles are determined analytically by the remaining k — 3 angles. 



2.2 Perturbations 

A polygon P m with m edges is determined up to rigid motion by m—1 consecutive edge 
lengths si, . . . , s m -i and the m— 2 consecutive interior angles fa, . . . , 4> m -2 with fa between 
edges Si and Sj+i. 

The space of perturbations of the angles of an m-gon P m which preserve the edge lengths 
is m— 3-dimensional, and is generated by "local" perturbations which change only four con- 
secutive angles. Here by perturbation we mean the derivative at of a smooth one-parameter 
path in the space of m-gons with the same edge lengths as P m . Such a perturbation is de- 
termined by the derivatives of the angles with respect to the parameter t along the path. 
We define to be the infinitesimal perturbation of the angles of P m , preserving the edge 

lengths, for which ^ = unless j is one of i — 1, i, i + 1, i + 2 (indices chosen cyclically) and 

f& = 1. See Figure 11 




Figure 4: Local perturbation of vertex 4 of an octagon 
For generic P, the -ir for % — 1, 2, . . . ,m — 3 generate all edge-length preserving perturba- 



at 



tions of P. These -M- are useful because they provide a local infinitesimal coordinate charts 
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for the boundaries of the various sets BP^(T) which share the same cycle P m . 

For example, suppose that BP R (T) for some T is parametrized by angles 9%, . . . , n _i, 
and we are on a part of the boundary defined by a cycle with interior angles <pi, . . . ,<p m 
(so the 0's are differences of the #'s). The infinitesimal (signed) volume of the part of the 
boundary swept out by the perturbations -J^-, for i = 1, . . . , m — 3 is 

(#x A • ■ ■ A rf0 m _ 3 ) ( • • • , = det . (1) 



That is to say, if P(ti, . . . , t m _ 3 ) for (ti, . . . , t m ^ 3 ) G [0, e]' m ~ 3 is a perturbation of P m , with 
P(0, • • • , 0) = P m and with d-P(o,.-.,o) having the above properties, then the signed volume 
swept out by all these perturbations, divided by e m ~ 3 , is given by (fH). 

We also need to consider perturbations of P m which change the edge lengths. Let ^ 
be the perturbation in which two consecutive edge lengths increase and the rest remain 
unchanged, and only three angles change: 

— % T = — ^- = 1 and — f = for other j, and 
obi obi 

— f = unless 7 = z — 1. i or z + 1. 

This is obtained by moving only vertex i. Note that ^ for i = l,...,n generate all 
infinitesimal perturbations of the edge lengths. With the above these generate all 
motions of P m . 



2.3 Volumes 

Here we determine how the volume of BPk(T) changes when one of the radii is increased. 

Let X be an order-n polymer in the boundary of BP^(T). We assume that X is in a 
codimension-1 part of the boundary, that is, it has only one cycle, C, with vertices v i, . . . , v m 
in counterclockwise order. Let <f>x, . . . , <p m be the corresponding interior angles of C. Let us 
assume that the tree T contains all edges of C except the edge between f m _i and v m , so that 
angles (pi, ... , 4> m -2 parametrize the polymers in BP^(T) close to X. 

Recall that each edge of T is directed away from the root. Some of these are directed 
clockwise around C and some counterclockwise. Let Vi be the vertex of C closest to the root. 
The edges of C are oriented clockwise around C to the left of Vi and counterclockwise to the 
right of Vi. Let 9 Vj be the angle of the edge of T whose head is at Vj. 

Lemma 1. 

d6 vi A • • • A A • • • A d9 Vm = (-l)™" 2 ^! A d<j) 2 A • • • A rf0 m _ 2 A d6 Vm . (2) 

Proof. We can rewrite the left-hand side in terms of the d(j)j as follows. To the right of 
Vi, we have <pj = n — (9 V . +1 — 9 V ). To the left of v i: we have (fik = it + Vk _ 1 — 6 Vk . Also, 
(pi = 9 Vi l — 9 Vi+1 . Replace successively d9 Vk by — d(p k for k = l,2,...,i — 1. Then replace 
successively d9 Ve by —d(fig-i for £ — i + 1, i + 2, . . . , m — 1. There are precisely m— 2 minus 
signs. □ 
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The volume form on BP^(T) is d9\ A ■ • • A d9 n -%. We can write this as a A (3, where a is 
the form on the left-hand side of ([2]). The form f3 involves the edges which are not part of 
the cycle C. 

If we fix the angles of a polymer X in the boundary of BP^(T) and change the radii by 
a small amount, from R to R', the boundary of BP#/(T) moves in general away from X. 
To compute the change in volume of BP^(T), we integrate, along the entire codimension-1 
boundary of BP^(T), this displacement times the volume form of the boundary. Let us con- 
sider a radius perturbation On the boundary we use the local infinitesimal coordinates 
defined by the perturbations -J^- for i = 1, . . . ,m — 3. The angles 9 which are not part of 
the cycle C can be perturbed independently; let J- be a perturbation of an angle 9j not 
on C, which changes only this angle. The angle 9 m can also be perturbed independently of 
the angles cj)j in the cycle and other angles: it just determines the orientation of the cycle. 
Let Tfjj— be a perturbation of 9 m . Let be the product of the perturbations for 9 m and the 
remaining angles not in C. 

The local volume element gained or lost by BP#(T) is then the product of /3(gg) and 

, f d d 9 \ 

u = {dcpi A ■ ■ • A d(f) m _2) ^r, • • • , ^- • 

\dSj dti dt m ^J 

The total volume change of BP^(T) is the integral of u times over the entire 

boundary of BP#(T). (More precisely, there is a corresponding form u = uj(C) for each 
codimension-1 piece of the boundary corresponding to the possible cycles C formed by T. 
The sum of the integrals of each form u>(C) times over the corresponding part of the 

boundary gives the total volume change.) 

Now as T ranges over the trees obtained by removing one edge from C, the individual 
forms on the right-hand side of ([2]) are obtained from d(f)\ A • • ■ A d<p m by removing two 
consecutive d&'s, and wedging the result with d9 Vm (which, since it plays the role of a global 
rotation of the cycle, is the same for each i). 

To prove that the sum of the volumes of the BP^(T) is constant, it suffices now to observe 
the following: 

Lemma 2. If 4>i + ■ ■ ■ + <j) m is constant, then 

2_. d(j)i A • • ■ A G?0j_i A d(j) i+ 2 A • • ■ A d(p m = 
i=i 

( with cyclic indices, and where if m is even we must put a — sign in front of the last term 
d<p 2 A • • - A d<f) m -i). 

Proof. Substitute d<p m = —difii — ■ ■ ■ — d<p m ^i in each term and simplify. □ 

It remains to show that this constant volume in fact takes the claimed value. 

Theorem 3. For any radius vector R of length n, the volume of the space of branched 
polymers is (n— l)!(27r)" _1 . 
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Proof. Choose e > very small and let R be given by r« = e l . Let X be a uniformly random 
configuration of disks with these radii, forming some labeled tree T. Suppose that for some 
j < n, disks 1 through j are connected. Then we claim that with probability near 1, disk j+1 
touches one of disks 1 through j. To see this, observe that otherwise disk j+1 is connected 
to some previous disk i, 1 < i < j, via a chain of (relatively) tiny disks whose indices all 
exceed j+1. Let disk k, k > j+1, be the one that touches disk then the angle of the 
vector from the center of disk k to the center of disk j+1 is constrained to a small range, 
else disk j+1 would overlap disk i. It follows that BP^(T) has lost almost an entire degree 
of freedom, thus has very small volume; in other words, the tree T is very unlikely. 

Suppose, on the other hand, that for every j, disks 1 through j are connected. Then we 
may think of X as having been built by adding touching disks in index order, and since each 
is tiny compared to all previous disks, there is almost a full range 2n of angles available to 
it without danger of overlap. 

It follows that as s — > the volume of the space of polymers with radius vector R 
approaches the volume of the space of crossing inductive trees, namely (n— l)!(27r) n_1 . Since 
this volume does not depend on R, we have equality. □ 

2.4 Generalization to graphs 

Let G be a graph on vertices {1, . . . , n} whose edges are equipped with positive real lengths 
Tij . A G-polymer is a configuration of points in the plane, also labeled by {1, . . . , n}, such 
that: 

1. point number 1 is at the origin; 

2. for each edge {i, j} of G, the distance p(i,j) between points % and j is at least r^; and 

3. the edges {i,j} for which p(i,j) = span G. 

We denote the set of G-polymers realizing a given (spanning) tree T by BPc(T). 

Note that if R = (r 1; . . . ,r n ), and G is the complete graph K n with = + rj, then 
a G-polymer is precisely the set of centers of the disks of a polymer with radius vector R, 
in the sense of the previous sections. The volume Vq of the space of G-polymers is defined 
as before by the angles made by the vectors from % to j, where {i,j} is an edge for which 

p(hf) = r ij- 

In fact, the proof of Lemmas [1] and [2] extend without modification to show that Vq does 
not depend on the lengths (even if they fail to satisfy the triangle inequality), but only 
on the structure of G. This leaves us with the question of computing Vq for a simple graph 
G. 

To do this, we label the edges of G arbitrarily as ei, . . . , e m and if = {i,j} we choose 
its edge-length to be e k for e > and very small. Then (since e < |), for the volume 
of BPq{T) to be non-zero, there must not be an edge of G \ T such that k is the lowest 
index of all edges in the cycle made by adjoining ek to T. If no such edge exists we say that 
T is "safe"; and in that case, arguing as in the proof of Theorem [31 there is almost no danger 
of violating condition (2) above in a random element of BPq(T). Thus the volume of the 
space of configurations in BPq(T) is nearly the full (27r) n_1 . 
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It follows that the volume of the space of all G-polymers is //(G)(27r)™ _1 , where //(G) is 
the number of safe spanning trees of G. Since u(G) does not depend on the edge labeling, 
one might suspect that it has a symmetric definition, and indeed it does. 

Lemma 4. For any graph G, the number //(G) of safe spanning trees of G is equal to the 
absolute value of the sum over all spanning subgraphs H of G, of (— ly H \. 

Proof. A simple inclusion-exclusion argument suffices. Let us fix any numbering of the edges 
of G and, for each spanning tree T, let B(T) be the set of "bad" edges of G \ T, that is, 
edges which boast the lowest index of any edge in the cycle formed with T. Associate to 
each spanning graph H the spanning tree T(H) obtained by repeatedly removing the lowest- 
indexed edge from each cycle. Then for n odd, a spanning tree T with set B(T) of bad edges 
is counted once positively for each even subset of B(T) and once negatively for each odd 
subset; and vice- versa for n even. It follows that in the sum (which we denote by //(G)) T 
has a net count of unless B(T) is empty, in which case it counts once positively (for n odd) 
or negatively (n even). But the trees for which B(T) is empty are exactly the safe trees. □ 

We conclude: 

Theorem 5. The volume of the space of G -polymers in the plane is fi(G)(27r) n ~ 1 . 

Comparing with Theorem [3], we have indirectly shown that n(K n ) = (— l) n ~ l (n— 1)!. In 
general //(G) = \Tg(0, 1)| where Tq is the Tutte polynomial of G (see e.g. [H HI [17]). We 
note also that //(G) plays the role of Brydges and Imbrie's function "Jc" in the dimension-2 
case. 

The computation of T G (0, 1), hence also of //(G), is unfortunately #P-hard for general G 
[9]. The point (0,1) is not, however, in the region of the plane in which Goldberg and Jerrum 
[6] have recently shown the Tutte polynomial to be hard even to approximate. Thus, there 
is some hope that a "fully polynomial randomized approximation scheme" can be found for 
//(G). 

We conclude this section with a new solution of a notoriously difficult puzzle, which 
appears as an exercise in p3], derived from Rayleigh's investigation (see [20]) of "random 
flight." The exercise calls for proving the corollary below by developing the Fourier analysis 
of spherically symmetric functions, then deriving a certain identity involving Bessel functions. 
Curiously, it is (we believe) the only mention of continuous random walk in Spitzer's entire 
book. 

Corollary 6. Let W be an n-step random walk in M. 2 , each step being an independent 
uniformly random unit vector. Then the probability that W ends within distance 1 of its 
starting point is 1/(71+1). 

Proof. The volume of the space of such walks, beginning from the origin, is of course (27r) n . 
If the walk does not terminate inside the unit disk at the origin, it is in effect a C n+ \- 
polymer, where G n +i is the n+1 cycle in which vertex % is adjacent to vertex z+1, modulo 
n+1. Since //(G n+ i) = |1 — (n+l)| = n, the volume of the space of G n+ i-polymers is n(27r) n . 
Since the spanning tree with no edge between nodes 1 and n+1 is one of n+1 symmetric 
choices, the volume of the G n+ i-polymers which correspond to non-returning random walks 
is n(27r) n /(n+l), and the result follows. □ 
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2.5 The bipartite case 



One special graph of interest is the complete bipartite graph K m ^ representing particles of 
two types, each particle interacting only with particles of the other type. Note that in the 
hard-core model, phase transition has been proved in this situation [13] — in dimensions 2 
and higher — but not for the complete graph. 

Put yU mn = fi(K m n ) and let H(x,y) be the exponential generating function for /i m>n , 
given by 

m=l n=l 

Lemma 7. 

H(x, y) = \og(e~ x + e~ y - e~ x - y ) . 

Proof. Observe that 

A*m,n = X^ _1 ) fc/i ( mi ' n i)M m 2, n 2 ) . . . t-i{m k , n k )C(mi, . . . , m k , n u . . . , n k ) 

where the sum is over all partitions m = m\ + • • • + m k and n—l = n\ + ■ ■ ■ + n k , and C is 
a combinatorial factor — the product of two multinomial coefficients divided by appropriate 
factorials when some of the pairs (m^, rii) are equal. 

We claim that dH(x,y)/dx = —1 + e~ y ~ H( - x ' y \ This is because dH/dx in effect removes 
one vertex from the left (m-vertex) side of G; this leaves singleton vertices from the right 
side, counted by e y , and other components, counted by e H ( x ' y \ The — signs are contributed 
by edges connecting these components to the removed vertex, and the —1 summand comes 
about because there must be at least one remaining component. 

Solving the differential equation with initial condition H (0, y) = yields 

H(x, y) = hg(e~ x + e~ y - e~ x ' y ) = -x - y + \og(e x + e y - 1) . 

□ 

The argument generalizes to the complete fc-partite case, giving 

k / k 

H(x u . . . , x k ) = - x i + lo S 1 - k + eXk 

i=l \ i=l 

2.6 Construction 

We now show inductively how to construct a uniformly random branched polymer of order 
n in the plane. 

We begin with a unit disk centered at the origin. Suppose we have constructed a polymer 
of size n—l, n > 1. We choose a uniformly random disk from among the n—l we have so 
far, then choose a uniformly random boundary point on that disk and start growing a new 
disk tangent to that point. If a disk of radius 1 fits at that point, this will define a polymer 
of size n. 
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Otherwise there is a radius < r < 1 at which a cycle forms with the new disk and 
some other disks present. At this point our polymer X is in the boundary of the space 
BP j r(T), where R — {1,1, . . . ,l,r}, and we need to choose some other tree T' for which X 
is in the boundary of BP^(T'), and which has the property that increasing r (and leaving 
the angles fixed) will not cause the disks to overlap. There will be at least one possible such 
T' because the volume of BP#(T) is decreasing as r increases and so must be compensated 
by an increase in volume of some BP^(T'). We choose randomly among the BP/j(T') with 
increasing volume, with probability proportional to the infinitesimal change in the volumes 
of the BP/j(T')'s as r increases. This ensures that the volume lost to BP#(T) as r increases 
is distributed among the other BP^(T') so as to maintain the uniform measure. (In the 
language of Markov chains, this is the detailed balance condition). 

Figure [5] shows snapshots of the construction of a random polymer, in the process of 
growing its third and fourth disks; Figure [6] shows a polymer of order 500 generated by this 
method. 

Figure 5: A random planar branched polymer growing new disks 

All of the above is easily generalized to produce uniformly random G-polymers for any 
connected graph G with specified edge-lengths (and in fact we will need this construction 
later, when generating 3-dimensional polymers). The vertices of G may be taken in any order 
v i, . . . , v n having the property that the subgraph Gk induced by Vi, . . . , is connected for all 
k. When a uniformly random Gfc_i-polymer has been constructed, a new point corresponding 
to vertex Vk is added coincident to a point uniformly chosen from its neighborhood — in other 
words, we start by assuming that the edges of Gk incident to v k are infinitesimal in length. 
These edges are then grown to their specified sizes, breaking cycles when they are formed in 
accordance with the rules above. 
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Figure 6: A uniformly random two-dimejnsional branched polymer of 500 disks. 



3 The 3-dimensional Case 



3.1 Volume invariance 

Branched polymers in 3-space share many of the features of planar branched polymers. [2] 
showed that the volume of the configuration space of polymers in 3-space is n n ~ l (27r) n ~ 1 . 
Whereas the planar configuration space volume was independent of the radii of the balls, the 
same is not true in 3 dimensions. However, there is an invariance inherited from the plane 
under a different change of norm, which we now describe. 

Let G be a graph with n vertices and edge weights /3y > 0. A 3-dimensional G-polymer 
is a set of n points v i, . . . , v n G M 3 such that for all i,j we have 

IK - vj f := (vj - v)f + frM ~ + K 3 - > h 

with equality holding on a spanning tree of G. 

When all (3^ are 1 this defines the standard branched polymer. Note that if ||t>j — 
Vj\\ = 1 then Vj is on the surface of an spheroid centered at V{. We measure the volume of 
the configuration space of 3-dimensional polymers using the normalized surface area of the 
corresponding spheroids; we will see that this volume is independent of the fy. 

3.2 One-dimensional projections 

Recall that the surface area measure of a sphere S 2 , projected to a line running through its 
center, projects to 2ir times Lebesgue measure on the image segment. The same is true of 
the spheroid {y G IR 3 : \\v\\ = 1} for any (3, when projected to the x-axis, and it follows 
that for purposes of computing the volume of the configuration space, we may assume that 
the polymers are parametrized by the length of the projection of each v% — Vj on the x-axis 
together with its angle to the positive y-axis when projected onto the yz-plane. 

Let xi, . . . , x n be the projections of V\, . . . , v n to the x-axis. We suppose, after relabeling 
if necessary, that the x, are ordered x\ < X2 < ■ ■ ■ < x n . If v i and Vj are adjacent in the 
polymer then \xi — Xj\ < 1. (See Figure 0) 

Lemma 8. Fix a graph G on {l,...,n} with edge weights {/%}■ The n—1- dimensional 
volume of the set of G -polymers whose centers project to x\ < ■ ■ ■ < x n is an integer multiple 
of (27r) n_1 and depends only on the set of pairs i,j with \xj — x.i\ > 1. 

Proof. In any such polymer, the distance between the yz-plane projections of each pair i,j 
of adjacent centers is some fixed depending only on \xi — Xj\ and /3^. For non-adjacent 
centers, this distance is at least some provided — Xj\ < 1; otherwise it is unconstrained. 

It follows that if we let H be the graph on vertices {1, . . . , n} given by i ~ j iff \x— Xj \ < 1, 
and define K to be the intersection of (the edges of) G and H (with edge-lengths supplied 
by G), then by Theorem [5] the desired volume is [i(K) ■ (27r)" _1 . □ 

3.3 Complete graph 

Note that if is a "unit interval graph" (see e.g. [12]) defined by overlapping unit-length 
intervals, in this case with their centers at the Xj. In the case of standard polymers, where 
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Figure 7: A branched polymer projected onto the x-axis and yz-p\ane. 

G is complete, K is just H with appropriate edge lengths assigned. Thus fi(K) = fi(H) and 
this is easy to compute: taking an arbitrary total order on the edges and /3 ek = e k for e small 
(where is the kth edge), the safe spanning trees of H are only those which are "inductive" 
in the sense of the Introduction: all paths from the root are increasing. It follows that each 
vertex j > 1 has as its parent some i < j for which X{ — Xj < 1, thus 

n 

ft 7 0') 

where is the number of % < j for which Xj — Xi < 1. 

It will be convenient temporarily to limit discussion to the space BP^ of order-n polymers 
for which X\ = 0, i.e. those whose roots extend farthest to the left along the x-axis. Define 
the type cr{X) of a polymer X in BP^ to be the permutation a of {2,3, . . . , n} for which 
x 'v{2) < x ct(3) < " " " < x ' a (n)i where x\:= Xi mod 1. Then fi(H) depends only on o and we 
may call it /x(cr). The projections uniquely determined by a and the arbitrary 

subset {x[, . . . , x' n } of [0, 1], thus: 

Corollary 9. The volume of the space of branched polymers of type o in BP n is precisely 

When a is the identity permutation /, all the Xi are in [0,1], and fx takes its maxi- 
mum value (n— 1)!. This is just the planar case in disguise (although note that the three- 
dimensional volume of the polymers of type / incurs another factor of n on account of the 
choice of k for which x^ = 0). On the other end of the scale, the minimum value fi(o~) = 1 
is achieved in the spread-out case when Xi + 2 — Xi > 1 for each i = 1, . . . , n— 2; the number 
of such a is the "Euler number" E n -\ (see e.g. Stanley [T5]). 
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Let T n be a uniformly random tree on the labels {1, . . . ,n}, with an independent uni- 
formly random real length U{j in [0,1] assigned to each edge For each j = 1, . . . , n let 
dj be the sum of the lengths of the edges in the path from the root (vertex 1) to j in T; 
and let = bi < b2 < ■ ■ ■ < b n be the dj taken in order. Let B be the (random) vector 
(6i, . . .A)- 

Theorem 10. Let X be a random branched polymer from BP n , and = X\ < X2 < • • • < x n 

the projections of its centers onto the x-axis. Then the random vector (0,a;2,X3, . . . ,x n ) is 
distributed as B. 

Proof. Suppose first that a tree T n is fixed and that its B-vector is of type a. If we allow the 
edge-lengths of T n to vary, we find that to maintain type a the quantities b<i mod 1, . . . , b n 
mod 1, which are independent, uniformly random drawings from [0,1], must fall in a partic- 
ular order. Thus the probability that the edge-length assignments to any particular combi- 
natorial tree T will yield a B-vector of any fixed type a is either or a constant independent 
of T and o. 

In view of Corollary [9J it suffices then to show that the number of labeled trees T n which 
contribute to type a is /i(cx), but this is easy. Given B, the node of T n corresponding to bj 
must have as its parent (counting node 1 as root) some node corresponding to an i < j for 
which bj — bi < 1. □ 

Theorem [H)] says that the x-axis projections of a random X e BP^ can be obtained by 
planting vertex 1 of T n at x = and stretching the tree to the right, letting the rest of its 
nodes mark the projections. 




6 6 3 1 



13 11 
12 10 
110 
10 10 
10 

10 



Figure 8: The matrix of types and trees for n — 4. 

Figure [8] illustrates the case n — 4. The rows are indexed by types, presented as sample 
projections, each accompied by its relative volume /x(cr). The columns are indexed by trees, 
each weighted by its number of distinct labelings (consistent with root at 1). 



4 ■ ■ I * 4 

— h — * 2 
<j> — • — 1 — • — • — 2 

<r— 1 ' — I 1 

* H — • — h* 1 



11 



Note that the theorem does not say that the tree structure of a random 3-dimensional 
polymer is uniformly random; for example, no polymer can have a node of degree greater 
than 12. It does imply, however, when combined with the proof of Theorem [8], that if the 
polymer is not made of spheres but of ellipsoids with widely ranging y — z axes, randomly 
assigned to labels, then indeed the tree structure approaches uniformly random labeled trees. 

From Theorem [10] we can incidentally deduce the not completely obvious fact that the 
"reverse" vector (0, b n — 6 n _i, b n — b n ^ 2 , ■ ■ ■ > b n } has the same distribution as B. For polymers, 
consequences of the theorem include the Brydges-Imbrie volume calculation and more: 

Theorem 11. The total volume of the space of 3-dimensional branched polymers of order 
n is n n ~ 1 (2ii) n ~ 1 , and the expected diameter (combinatorial or Euclidean) of a random such 
polymer grows as n 1 ^ 2 . 

Proof. For the volume, we apply Cayley's theorem (to the effect that the number of labeled 
n-node trees is n n ~ 2 ) and the fact that, since relabeling of polymers preserves volume, the 
volume of BP^ is just 1/n times the volume of the whole space. 

For the diameter we make use of Szekeres' Theorem (see [HI [16]) saying that the expected 
length of the longest path in a random tree on n labels is of order y/n. The expected length 
of the longest path from the root in our edge-weighted tree T n must therefore also be of 
order y/n, and this is exactly the length of the projection of our random polymer on the 
a;-axis. Since the space of polymers is independent of choice of axes, the spatial diameter of 
a random polymer must also be of order \fn. □ 

3.4 Construction 

To construct a uniformly random three-dimensional branched polymer of order n, we first 
select a uniformly random labeled tree T, then a set {x±, . . . , x n } of projected centers on the 
x-axis. We then build the yz-plane projection using our 2-dimensional polymer construction; 
this yields the locations of the n centers in 3-space, and it remains only to pick a root and 
translate it to the origin. 

The tree T on vertices {1, . . . , n} can be selected from the n n ~ 2 possibilities by means of 
a Priifer code (see, e.g., [8]), which is itself just a sequence of n— 2 numbers between 1 and 
n. The first entry of the code is the label of the vertex adjacent to the least-labeled leaf of 
T; that leaf is then deleted and succeeding entries defined similarly. The reverse process is 
also unique and easy. 

The projections are defined by assigning independent uniformly random reals u G [0, 1] to 
each edge of T, then letting X{ be the length of the path from vertex i to vertex 1. The unit- 
interval graph H is defined as above on the tree- vertices, namely by % ~ j if \xj — Xi\ < 1. 
Edge-lengths are assigned to H by t{i, j) = a/1 — (x~j — Xi) 2 so that the spheres of the 
polymer corresponding to tree vertices i and j are touching just when their centers lie at 
distance £(i,j) when projected onto the T/z-plane, and in any case lie at least that far apart. 

From the argument above we know that given x\, . . . ,x n , the yz-pl&ne projections are 
exactly a uniformly random planar if -polymer, which we then select using the methods of 
Section 12.61 

Combining the x-axis and ?/z-plane projections gives us the centers of a uniformly random 
branched polymer in 3-space (with spheres of diameter 1), except that sphere number 1 is 
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forced to have its center on the yz-p\ane; we now choose a sphere uniformly at random to 
be the new root, and translate the polymer so that this sphere's center is at the origin. 

Figures M, HQl and ?? are snapshots, from three angles, of a 3-dimensional branched 
polymer constructed as above. 

4 Open problems 

1. Is there a geometric interpretation of the local volume changes of the BP^(T)-which 
clearly depends on the shape of the cycle C? This would lead to a possible natural 
geometrization of the space of polymers. 

2. What are the volumes of BP#(T) for each T? 

3. What is the expected diameter (combinatorial or geometric) of a random two-dimensional 
branched polymer? 

4. More generally, what do random polymers look like in the scaling limit, in any fixed 
dimension? 
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Figure 10: The same polymer, slightly rotated 
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